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Abstract 

Context. Recently X-ray emission from protostellar jets has been detected with both XMM-Newton and Chandra satellites, but the physical 
mechanism which can give rise to this emission is still unclear. 

Aims. We performed an extensive exploration of a wide space of the main parameters influencing the jet/ambient interaction. Aims include: 1) 
to constrain the jet/ambient interaction regimes leading to the X-ray emission observed in Herbig-Haro objects in terms of the emission by a 
shock forming at the interaction front between a continuous supersonic jet and the surrounding medium; 2) to derive detailed predictions to be 
compared with optical and X-ray observations of protostellar jets; 3) to get insight into the protostellar jet's physical conditions. 
Methods. We performed a set of bidimensional hydrodynamic numerical simulations, in cylindrical coordinates, modeling supersonic jets 
ramming a uniform ambient medium. The model takes into account the most relevant physical effects, namely the thermal conduction and the 
radiative losses. 

Results. Our model explains the observed X-ray emission from protostellar jets in a natural way. In particular we find that the case of a 
protostellar jet less dense than the ambient medium reproduces well the observations of the nearest Herbig-Haro object, HH 154, and allows us 
to make detailed predictions of a possible X-ray source proper motion (v B h ~ 500 km s _1 ), detectable with Chandra. Furthermore our results 
suggest that the simulated protostellar jets which best reproduce the X-rays observations cannot drive molecular outflows. 
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1. Introduction 

The early stages of the star birth are characterized by a variety 
of mass ejection phenomena including collimated jets. These 
plasma jets can travel through the interstellar medium at super- 
sonic speed with shock fronts forming at the interaction front 
between the jet and the unperturbed ambient medium. In the 
last 50 years these features have been studied in detail in the 
radio, infrared, optical and U V bands and are known as Herbig- 
Haro (hereafter HH) obj ects dHerbigll950HHaroll952t see also 
iReipurth & BalhfcOOll) . 

Pravdoetal predicted that the most energetic HH 

objects could be sources of strong X-ray emission. Following 
IZeFdovich & Raizerl i 19661) one can derive useful relations be- 
tween the physical parameters of interest (as the plasma tem- 



Send offprint requests to: R. Bonito 
e-mail: sbonito@astropa.unipa.it 



perature and the shock velocity) of the post-shock region, in 
particular 

T - 7-1 f mV sh\ m 

Tpsh - w+w V^J (1) 

where T ps h is the post-shock temperature, 7 is the ratio of spe- 
cific heats, t? s h is the shock front speed, m is the mean parti- 
cle mass and k-g, is the Boltzmann constant. Assuming a typi- 
cal velocity, v*h ~ 500 km s _1 (as measured in HH 154, see 
iFridlund et al.l2005l) . the expected post-shock temperature is of 
few millions degrees, thus leading to X-ray emission. 

Recently, X-ray emission from HH objects has been de- 
tected with both the XMM-Newton and Chandra satellites: 
the low mass you ng stellar objects (YSO) HH 2 in Orion 
dPravdo etaj 2001) and HH 154 in Taurus JFavata et all2 002: 
lBallvetal.1120031). the high ma ss YSO objects HH 80/81 in 
Sagittarius (Prav doetal J 12004^ and HH168 in Cepheus A 
dPravdo & Tsuboill2005h . and HH 210 in Orion dGrosso et alJ 
2006j). Indications o f X-ray emission from prot ostellar jets 
are also discussed by iTsuiimoto et alJ d2004l) and iGiidel et all 
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Table 1. Relevant physical quantities observed in confirmed X-ray emitting HH objects, where Lx is the reported X-ray lumi- 
nosity, kT and Nr are the best fit parameters derived from spectral analysis for the temperature (in keV) and for the interstellar 
absorption column density, is the shock front velocity derived from opti cal observations, D is the distance of the object 
observed and Lh^/Lp, is the bolometric luminosity of the jet drivin g source ( Pravdo et alJ2 001 : Fav ata et alJ l2002: Ball v et alJ 
l2003UPravdo et alJ2004trPravdo & TsuboJ2005tlcir"osso et alJ200l . 



object 


L x 


kT 


N H 




D 


Lhol/ Lq 


Lx/L bo i 




[10 29 ergs" 1 ] 


[keV] 


[10 22 cm" 2 ] 


[km s" 1 ] 


[pc] 






HH 2 


5.2 


0.23 


< 0.09 


230 


480 


8r 


1.7 x 10" b 


HH154 


3.0 


0.34 


1.40 


500 


140 


40 b 


9.7 x 10 -7 


HH 80/81 


450 


0.13 


0.44 


700 


1700 


2 x 10 4c 


1.5 x 10" 4 


HH 168 


1.1 


0.5 


0.40 


500 


730 


2.5 x 10 4c 


3.6 x 10~ 7 


HH210 


10 


0.07 - 0.33 


0.80 


130 


450 







a Chini et al. (2001) 
b Liseau et al. (2005) 
c Curiel et al. (2006) 



J2005h . A summary of the relevant physical quantities observed 
for these objects is presented in Tab.^ 

In addition to the intrinsic interest in their physics, under- 
standing the X-ray emission from protostellar jets is important 
in the context of the physics of stars and planets formation. 
X-rays (and more in general ionizing radiation) affect many 
aspects of the environment of young stellar objects and, in par- 
ticular, the physics and chemistry of the accretion disk and its 
planet-forming environment. The ionization state of the accre- 
tion disk around young stellar objects will determine its cou- 
pling to the ambient and protostellar magnetic field, and thus, 
for example, influence its turbulent transport. In turn, this will 
affect the accretion rate and the formation of structures in the 
disk and, therefore, the formation of planets. Also, X-rays can 
act as catalysts of chemical reactions in the disk's ice and dust 
grains, significantly affecting its chemistry and mineralogy. 

The ability of the forming star to ionize its environment will 
therefore significantly affect the outcome of the process, inde- 
pendently from the origin of the ionizing radiation. While all 
young stellar objects are strong X-ray sources, they will irra- 
diate the disk from the central hole, so that stellar X-rays will 
illuminate the disk with grazing incidence, concentrating their 
effects in the central region of the disk (although flared disks 
can alleviate the problem to some extent). Protostellar jets are 
located above the disk, so that they will illuminate the disk 
with near normal incidence, maximizing their effects even in 
the outer disk regions shielded from the stellar X-rays. For ex- 
ample, the emission from HH 154 is located at some 150 AU 
from the protostar, ensuring illumination of the disk with favor- 
able geometry out to few hundreds AU. 

Several models have been proposed to explain the X-ray 
emission from proto stellar jets, but the actual emission mecha- 
nism is still unclear. iBallv et all J2003) speculated on different 
mechanisms for the X-ray emission from HH 154: X-ray emis- 
sion from the central star reflected by a dense medium, X-ray 
emission produced when the stellar wind shocks against the 
win d from the com panion star, or produced in shocks in the 
jet. lRaga et al.N2002t) derived a simple analytic model, predict- 
ing X-ray emission originating from protostellar jets with the 
observed characteristics. 



Prompted by this recent detection of X-ray emission from 
HH objects, we developed a detailed hydrodynamic model of 
the interaction between a supersonic protostellar jet and the 
ambient medium, to explain the mechanism causing the X- 
ray emission observed. Our model takes into account opti- 
cally thin radiative lo sses an d thermal conduction effects. We 
use the FLASH code dFrvxell et al.l200oT) with customized nu- 
merical modules th at treat optically thin radiative losses and 
thermal conduction (Orland o et all2005l) . The core of FLASH 
is based on a directionally split Piecewise Parabolic Method 
(PP M) solver to handle c ompressible flows with shocks 
(IColella & Woodward fl984l) . FLASH uses t he PARAMESH 
librar y to handle adaptive mesh refinement (Mac Neice et alJ 
2000) and the Message-Passing Interface library to achieve par- 
allelization. 

In a previous paper jBonito et alJl2004T) . we presented a 
first set of results concerning a jet less dense than the ambient 
medium, with density contrast v = n a /rij = 10 (where n a is 
the ambient density and nj is the density of the jet) which emits 
X-rays in good agreement w i th the X-ray e mission observed in 
HH 154 jFavata et al.ll2002l) . bonito et alJ d2004 have shown 
the validity of the physical principle on which our model is 
based: a supersonic jet traveling through the ambient medium 
produce a shock at the jet/ambient interaction front leading to 
X-ray emission in good agreement with observations. In the 
present paper, we study the effects on the jet dynamics of vary- 
ing the parameters, such as the ambient-to-jet density ratio, 
v = n a /rij, and the Mach number, M = v- 3 /c a , through a 
wide range, to determine the range of parameters which can 
give rise to X-ray emission consistent with observations. Note 
that we use this definition for the Mach number to be able to 
compare the jet velocity with the ambient sound speed, to have 
information on how much the jet is supersonic. 

The paper is structured as follow: Sect. |2] describes the 
model and the numerical setup; in Sect.|5]we discuss the results 
of our numerical simulations; finally Sect.|4]is devoted to sum- 
mary and conclusions. In Appendix |X]we discuss our method 
to synthesize X-ray emission from our numerical simulations. 
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2. The model 

We model the propagation of a constantly driven protostellar jet 
through an isothermal and homogeneous medium. We assume 
that the fluid is fully ionized and that it can be regarded as a 
perfect gas with a ratio of specific heats 7 = 5/3. Also we 
assume a negligible magnetic field. 

The jet evolution is described by the fluid equations of 
mass, momentum and energy conservation, taking into account 
the effects of radiative losses and thermal conduction 



dp 
di 



+ V • p\ = 



+ V • p\\ + Vp = 

at 



dpE 



V • (pE +p)\ = -V -q- n e n H P(T) 



(2) 



(3) 



(4) 




Figure 1. Initial jet velocity (continuous line) and density 
(dashed line) as a function of the distance from the axis, r, 
for the light jet case discussed in Sect. 13.21 with parameters: 
M = 300, v = 10, rij 
w — 4 (see Tab. [5}. 



= 500 cm~ 6 , T a = 10 d K, r,- = 30 AU, 



where t is the time, p the mass density, v the plasma veloc- 
ity, p the pressure, q the heat flux, n e and njj are the elec- 
tron and hydrogen density respectively, P(T) is the optically 
thin radiative losses function per unit emission measure (for the 
P(T) we use a functional form, which takes into account: free- 
free, bound-free, bound-bound and 2 photons emission, see 
Raymond & Smithll 19771 iMewe et aill 19851 iKaastra & Mewel 
2000), T the plasma temperature, and 



(5) 



where E is the total energy and e the specific internal energy. 
We use the equation of state for an ideal gas 



P = (7 - 



(6) 



Following lPalton & Balbul i 19931) . we use an interpolation ex- 
pression for the thermal conductive flux of the form 



1 

9spi 



1 

9sat 



(7) 



which allows for a smooth transition between the classical and 
saturated conduction regime. In the above e xpression, g sp i rep- 
resents the classical conductive flux (Spitzer 1962|) 



-k(T)VT 



where k(T) = 9.2 x 10" 7 T 5 / 2 erg s" 1 K" 1 cm" 1 is the ther- 
mal c onductivity. The saturated flux, g sat , is JCowie & McKed 
1 19771) 



9sat 



-sign(S7T)5(f)pCs 



(9) 



where ~ 0.3 JGiulianll984tlBorkowski et all 19891 and ref- 
erences therein) and c s is the isothermal sound speed. 



2.1. Numerical setup 

We adopt a 2-D cylindrical (r, z) coordinate system with the 
jet axis coincident with the z-axis. For the different cases ana- 
lyzed, we have chosen different ranges for the radial and lon- 
gitudinal dimensions of the computational grid to follow in all 
cases the jet/ambient interaction for at least 20-50 years: the 
computational grid size varies from w 300 AU to w 600 AU in 
the r direction and from w 6000 AU to « 3 x 10 4 AU in the z 
direction. 

In the case of the jet less dense than the ambient medium 
(hereafter "light jet") which best reproduces observations, the 
integration domain extends over 300 AU in the radial direc- 
tion and over 6000 AU in the z direction. In the case of the 
jet with the same initial density as the ambient medium (here- 
after "equal-density jet") which best reproduces observations, 
the domain is (r X z) w (600 x 6000) AU. In this case, the 
radial axis is twice as large than in the light jet case because 
the cocoon surrounding the equal density jet has a radial exten- 
sion greater than in the light jet. The dimension of the compu- 
tational domain in the case of the jet denser than the ambient 
(hereafter "heavy jet") which best reproduces observations is 
(r x z) w (700 x 27000) AU. 

In all the cases, the initial jet velocity is along the z axis, 
coincident with the jet axis, and has a radial profile of the form 



V(r) 



V 



v cosh (r/r j ) w — (y — 1) 



(10) 



(8) where Vq is the on-axis velocity, v is the ambient to jet density 
ratio, r\ is the initial jet radius and w = 4 is the steepness 



parameter for the shear layer (as an example, see continuous 
line in Fig.^ f° r the light jet case discussed in Sect. l3.2l with 
parameters in Tab.|3j, to have a smooth transition of the kinetic 
energy at the interface between the jet and the ambient medium. 

The density variation in the radial direction (dashed line in 
Fig -ID is 



p(r) = Pj [v 



cosh(r/r,) 1 



(11) 
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where p\ is the jet density dBodo etalJl994 . 

Reflection boundary conditions are imposed along the jet 
axis, inflow boundary conditions at z = and r < rj and 
outflow boundary conditions elsewhere. 

The maximum spatial resolution achieved in the best light 
jet case, (in both r and z directions), is w 1.3 AU according to 
the PARAMESH methodology, using 4 refinement levels, cor- 
responding to covering the jet radius with 25 points at the max- 
imum resolution. The spatial resolution achieved in the equal- 
density case is half the one in the light jet case. In the best heavy 
jet model, the spatial resolution achieved is 8 times lower than 
in the light jet case. 

Our choice of different spatial resolution for the three cases, 
aimed at reducing computational cost, was necessary because 
the thermal conduction is solved explicitly in FLASH and, 
therefore, a time-step limiter depending on density, p, tempera- 
ture, T, and spatial resolution, Ax, is required to avo id numer- 
ical instability (see, for instance. lOrlando et all2 005). Stability 
is guaranteed for At < 0.5 Ax 2 /D, where D is the diffusion 
coefficient, related to the conductivity, k, and to the specific 
heat at constant volume, c v , by D = n(T)/(pc v ). In models 
characterized by high values of temperature (as, for instance, 
in the heavy jet case), therefore, a lower spatial resolution was 
required to avoid a very small time-step, At. 

2.2. Time scales 

Condensations of plasma, due to radiative cooling effects, can 
become thermally unstable; however the presence of thermal 
conduction can prevent such instabilities. By comparing the 
radiative, r ra d, and thermal conduction, r con d, characteristic 
times 



Trad 



Tcond 



P 



( 7 - l)n e n H P(T) 

p 7 I 2 _ 
% - 1 2 k(T)T ~ 



2.5 x 10 3 



1.5 x 10" 



T 3/2 



T 5/2 



(12) 



(13) 



where / represents the characteristic length of temperature vari- 
ations, we can infer which of the two competing processes 
dominates during the jet/ambient interaction. From the condi- 
tion 



( T ra d 
\ Tcond 



1/2 



(14) 



we ca n derive the cutoff length scale for instability, l-p dFiel 
1196.4 . which indicates the maximum length 



T 2 

l F » 1.3 x 10 6 — 

n 



(15) 



over which thermal conduction dominates over radiative ef- 
fects, in the classical conduction regime. An analogous estimate 
in the saturation regime leads to 



(M . « 3 x 10 7 



n 



which is one order of magnitude longer than the characteristic 
length in the classical regime. As discussed later in Sect. 13.21 
the comparison between the classical Field length (the short- 
est characteristic length) and the size of the region behind the 
shock at the head of the jet will allow us to determine if this 
region is thermally stable or not. 

In order to verify our hypothesis of a fully ionized gas, 
we computed the ionization time scale of the most relevant 
elements in the X-ray spectrum of a shocked plasma at T = 
3.4 x 10 6 K, assuming a post-shock density of about 10 4 cm -3 
(the light jet case). As an example, we derived that the ion- 
ization time scale for C and O is 1 to 2 orders of magnitudes 
smaller than the radiative and thermal conduction time scales 
so that the plasma can be considered in equilibrium. 

2.3. Parameters 

Our model solutions depend upon several physical parame- 
ters as, for instance, the jet and ambient temperature and den- 
sity, the jet velocity and its radius. With the aim to reduce the 
number of free parameters in our exploration of the parame- 
ter space, we fixed the jet radius to rj « 30 AU, according to 
iFavata et alJ J2002I) who found this characterist ic linear scale, 
from t he X-ray thermal fit, and according to Fridlund et al. 
(2005) who showed HST images of the internal knots of HH 
154 with dimension r w 30 AU at the base of the jet 1 . However 
detailed simulations with different rj values are not necessary 
since we can predict the effects of varying the jet radius from 
the model results we obtained so far. In fact we expect the X- 
ray emitting region to grow in size as rj grows. Since the X-ray 
luminosity is defined as Lx = n 2 VP(T), it depends on the 
cube of the radius. This means that, Lx being constrained from 
observations, a jet with a greater radius needs a lower density in 
order to reproduce observations. We impose an initial jet length 
Zj = 300 AU to avoid the ejected plasma that travels back in- 
side the boundary during the jet evolution. This choice of a 
non-zero initial jet length allows us to obtain an unperturbed 
boundary surface at z — 0. In all our simulations, we model a 
jet with initial density and temperature nj = 500 cm ~ 3 and 
Tj = 10 4 K r espectively, according to the values derived from 
observations (Fridlund & Liseaulll998l and lFavata et alJ l2002). 
The density and temperature of the ambient medium, n a and T a 
respectively, are derived from the choice of the ambient-to-jet 
density contrast, v and from the hypothesis of initial pressure 
balance between the ambient medium and the jet. We are left, 
therefore, with two non-dimensional control parameters: the jet 
Mach number, M, and the ambient-to-jet density contrast, v. 
For a more extended exploration of the parameter space, see 
Sect. 13.51 concerning the variation of the initial jet density, n- 3 . 
In our simulations we account for a wide jet/ambient parame- 
ters range shown in Tab.|2] 



(16) 



1 In Fridlund et al. (2005), on page 993 the authors discuss the 
working surface. The radius quoted is the one of the elongated Mach 
disk (probably representative of the jet) and it is 0.3/2" pb 30 AU. 
The separation between Mach disk and working surface is 0.6" or 4 
times this ~ 100 AU (M. Fridlund, private communication). 
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Table 2. Ra nge of parameters used in our numerical model (column 2) compared with typical outflow parameters (column 3 and 
4) shown in Bally & Reipurthl faOOA Tab. 1 . v is the ambient to jet density contrast; M is the Mach number; T a and n a are the 
ambient temperature and density, respectively; V\ is the initial jet velocity; v s h is the shock velocity of the knots inside the jet; M 
is the mass loss rate; L mcc h is the mechanical luminosity. 



Parameter Model Low Mass™ High Mass™ Units 

~v 0.01 ™ 300 

M 1^1000 

T a 30 ^10 6 - - K 

n a 5 10 5 - - cm -3 

Vj 85-^ 8500 - - kms" 1 

v sh 100- 2000 100^300 100 1000 kms" 1 

M 10~ 10 -J- 10" 8 10" 9 ~r~ 10~ 5 10" 6 ~r~ 10~ 2 M yr" 1 

imcch 6.7 x 10~ 5 ^ 67 0.001 -j- 1 0.1 -5- 1000 L Q 



a Bally & Reipurth (2002) 



In Sect. |5] we discuss the results derived from the explo- 
ration of the parameters space defined by M and v. 

3. Results 

3. 1. Exploration of the parameter space 

We performed a wide exploration of the parameter space de- 
fined by two free parameters: the jet Mach number, M — v-J c a , 
and the ambient-to-jet density ratio, v = n a /n- 3 (see Sect. l2.3> . 
The aim is to determine the range of parameters leading to X- 
ray emission from protostellar jets in agreement with the ob- 
servations. 

We first analyzed adiabatic hydrodynamic models, i.e. 
without thermal conduction and radiative losses. Then, for the 
most promising cases (i.e. those which reproduce the values 
of jet velocity, temperature and luminosity of the X-ray source 
derived from the observations), we performed more realistic 
simulations in which we have taken into account thermal con- 
duction and radiative losses effects. By comparing these mod- 
els with those without thermal conduction and radiative cool- 
ing, we explored how the presence of these physical processes 
affects the jet/ambient system evolution. We found that, in 
general, models with thermal conduction and radiation reach 
lower temperatures (up to 5 times lower than those achieved 
in the adiabatic cases). We also found that thermal conduction 
smooths the structures (well visible in the pure hydrodynamic 
cases) in the density and temperature distributions. 

In the following subsections, we discuss the models (shown 
in Fig. |2j in which both radiative losses and thermal conduc- 
tion are taken into account. In Fig.[2]green and red dots refer to 
those cases with X-ray luminosity Lx > 10 28 erg s~ x , shock 
front velocity v s h > 100 km s _1 and fitting temperature T < 
10 7 K, consistent with observations. We have chosen Lx one 
order of magnitude lower than the minimum value observed 
(see Tab. ^ to take into account fainter sources not detected 
so far; the r ed dot refers t o the representative case of HH 154 
discussed in lBonito et alJ J20 04). Squares show cases with ve- 
locity in the range of values observed, but with Lx < 10 28 erg 
s~ x . Diamonds mark the cases with velocity and X-ray lumi- 
nosity not consistent with observations. Triangles mark cases 
with temperatures higher than 10 7 K. The lower panel of Fig. 



Table 3. Summary of the initial physical parameters character- 
izing the "best-fit" models in the cases of light, equal-density 
and heavy jets: ambient-to-jet density contrast, v, jet Mach 
number, M, initial jet velocity, Vj, ambient density and the tem- 
perature, n a , and T a , respectively. In all the models, the initial 
jet density and temperature are rij = 500 cm ~ 3 and 2] = 10 4 
K, respectively. 



model 


V 


M 


«j 


n a 


Ta 








[kms^ 1 ] 


[cm~ 3 ] 


[10 4 K] 


light 


10 


300 


1400 


5000 


0.1 


equal-density 


1 


100 


1500 


500 


1 


heavy 


0.03 


30 


2500 


17 


30 



|2] show the initial velocity assumed in our simulations vs. the 
density contrast. 

From our exploration of the parameter space, we derived 
that the models in agreement with observations are included in 
a well constrained region. In the following sections, we discuss 
in details the "best-fit" models, i. e. those models which repro- 
duce X-ray luminosity and shock front speed values as close as 
possible to those observed, in the cases of light, equal-density 
and heavy jets (see Tab.|5J. 

3.2. Hydrodynamic evolution 

In Fig. [3] we show the mass density and temperature distribu- 
tions 20 years since the beginning of the jet/ambient medium 
interaction for the three best-fit models in Tab. [3] The light jet 
case is the one which better re produce the physica l param eters 
deri ved from obse rvations by Fridlu nd & Liseaul <|l 998) and 
Favat a et al for the HH 154 protostellar jet; its prop- 

erties have been discussed in lBonito et alJ (E004). 

In all the cases, at the head of the jet there is clear evidence 
of a shock front due to the plasma propagating supersonically 
along the jet axis. Just behind the shock front there is a local- 
ized hot and dense blob (see, for instance, the enlargement in 
Fig-EJfoi" the light jet case). 

The light jet is enveloped by a cocoon with temperature 
T w 7 x 10 5 K, almost uniform due to the thermal conduc- 
tion diffusive effect; nevertheless the cocoon temperature is 
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Figure 2. Exploration of the parameter space: jet Mach num- 
ber, M (upper panel) and initial jet velocity (lower panel) as 
a function of the ambient-to-jet density contrast, v. Green and 
red dots refer to models consistent with observations for X-ray 
luminosity, shock velocity and fitting temperature values (see 
text); diamonds refer to models which cannot reproduce obser- 
vations; squares refer to models in good agreement with obser- 
vations for shock velocity but not for X-ray luminosity values; 
triangles refer to models with too high (T > 10 7 K, one order 
of magnitude higher than observed) temperature. Empty sym- 
bols refer to pure hydrodynamic simulations; filled symbols re- 
fer to models in which thermal conduction and radiative losses 
effects are taken into account. 



not constant in time but decreases as the evolution goes on, 
leading to the formation of a cool and dense external enve- 
lope. Fig. |5] shows cuts of the density (continuous line) and 
temperature (dashed line) along the radius at z w 5000 AU, 
corresponding to the blob position 40 years since the begin- 
ning of the jet/ambient interaction: the hot (few millions de- 
grees) and dense blob is evident for r < 10 AU. The density 
decreases moving away from the jet axis along the radial direc- 
tion and, then, it increases again at the position corresponding 
to the external part of the cocoon. On the other hand, the tem- 
perature monotonously decreases, moving away from the jet 
axis along the radial direction. The blob, therefore, is expected 
to be an X-ray source; in Sect. 13.4.21 we will show that the 
X-ray source has luminosity and spectral characteristics con- 
sistent with those observed. 
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Figure 3. Mass density (upper semi-panels) and temperature 
(lower semi-panels) bi-dimensional cuts in the r — z plane after 
20 years since the beginning of the jet/ambient interaction for 
the best cases of light (upper panels), equal-density (middle 
panels) and heavy jet (lower panels). 
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Figure 4. An enlargement of the post shock region in the light 
jet case of temperature (left) and density (right) in linear scale, 
about 20 years since the beginning of the jet/ambient interac- 
tion. Note the hot and dense region just behind the shock front. 



The central panel in Fig.|5]shows 2-D sections in the (r, z) 
plane of the mass density and temperature distributions for 
the best-fit equal-density model (see Tab. The interaction 
between the protostellar jet and the ambient medium causes 
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Figure 5. Density (continuous line) and temperature (dashed 
line) profile along the radial direction at z w 5000 AU corre- 
sponding to the blob position at t = 40 yr (M = 300, v = 10). 

the presence of a dense and hot cocoon (n w 1400 cm -3 ; 
T w 2 x 10 6 K) surrounding the jet. Once again the cocoon 
is almost uniform for the presence of the thermal conduction 
but its temperature decreases with time. The cocoon becomes 
gradually cool and dense with time as in the light-jet case. Also 
in this case, the post-shock region is a hot and dense blob from 
which the X-ray emission originates (see Sect. l3.4.H for more 
details). 

In the heavy -jet case (lower panel in Fig. [3}, the jet is sur- 
rounded by a cocoon well smoothed by the effects of the ther- 
mal conduction and with a radial extension larger than in the 
other two cases. The cocoon has temperature of a few millions 
degrees and density lower than that of the jet. 

For the three best-fit models discussed above, we analyzed 
the thermal stability of the hot and dense blob localized be- 
hind the shock by comparing the size of the blob with the 
Field length Eq. ^] In the light jet case, the values obtained 
for the average blob temperature, T rj 3.4 x 10 6 K, and den- 
sity, n « 6500 cm~ 3 , lead to ~ 100 AU. Since the blob size, 
almost equal to twice the initial jet radius r- 3 w 30 AU (see Fig. 
|4}, is smaller than lp, it turns out that the blob is thermally 
stable. In the equal-density jet case, the density and tempera- 
ture of the blob at the head of the jet are n k, 1700 cm~ 3 and 
T » 4.3 x 10 6 , respectively, leading to l F « 10 3 AU. Also 
in this case, therefore, the blob is thermally stable, being its 
size ~ 100 AU, i.e. 10 times smaller than the Field length. In 
the heavy-jet case, the temperature and density of the blob are 
T ~ 10 7 K and n ~ 10 2 cm" 3 , leading to l F « 10 4 AU. Since 
the blob behind the shock front extends over about 100 AU, 
also in this case it is thermally stable. 

The position of the shock front as a function of time for 
the three best-fit models in Tab. [3] is shown in Fig. [6] For the 
light jet case, we derived an average shock velocity u s h « 500 
km s -1 , about 3 times lower than the initial jet velocity. This 
shock velocity is in good agreement with observed speeds in 
HH objects and in particular with that derived from HH 154 
data. Taking into accoun t the jet inclination w 45 degrees 
jFridlund & Liseaull 19981) . v s h ~ 500 km s~ 4 corresponds to 
a proper motion of » 350 km s — 1 , which, at the distance of 
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Figure 6. Shock front position vs. time. Dots mark the light jet 
case, stars the equal-density jet and crosses the heavy jet case. 
The lines refer to the best-fit from which we derive the average 
shock speed: v s h « 500 km s _1 , « 600 km s _1 , ks 1900 km 
s _1 for the light, equal-density and heavy jet case, respectively. 



HH 154, can be measured with well time-spaced Chandra ob- 
servations. 

In the equal-density jet scenario, we deduced an aver- 
age shock velocity v s h ~ 600 km s -1 , slightly greater 
than the value observed in HH 154 ( Fridl und & Liseaulll998t 
Favata et al I2002I) . but consistent with values observed in other 
HH objects (see Tab. [Q. For the heavy jet case, the average 
value of the shock speed is u s h ~ 1900 km s _1 which is too 
high with respect to the HH shock front velocities observed (cf. 
Tab.G). 

3.3. Emission measure distribution vs. temperature 

We derived the distribution of emission measure vs. tempera- 
ture, EM(T), in the temperature range [10 3 - 10 8 ] K at dif- 
ferent stages of the evolution of the jet/ambient system (see 
Appendix lAl for more details). Fig. shows the EM(T) for 
the three best-fit models in Tab. [3] after 20 years since the be- 
ginning of the jet/ambient interaction. 

In all the cases, we found that the shape of the EM(T) 
is characterized by two bumps, and does not change signifi- 
cantly during the system evolution. The relative weight of the 
bumps is different in the three cases. In the light jet case (upper 
panel in Fig. 0, the bumps are quite broad, the first centered 
at temperature T ~ 10 4 K with EM ~ 10 55 cm~ 3 , and the 
second one centered at T ~ 10 6 K with EM ~ 10 52 cm~ 3 , 
about three orders of magnitude lower than the previous one; 
the EM decreases rapidly above few millions degrees. 

In the equal-density jet case (middle panel in Fig. Q, the 
first bump is centered at T ~ 10 4 K with EM > 10 54 cm~ 3 , 
whereas the second bump is centered at T ~ 10 6 K with 
EM ~ 10 52 cm -3 , as in the light jet case. In the heavy jet 
case (lower panel in Fig. [7}, the EM(T) distribution appears 
flat with two weak peaks: the first centered at T ~ 10 4 K and 
the second one at a few millions degrees. Note that, in the heavy 
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Figure 7. Emission measure, EM, as a function of the temper- 
ature, T for the three models discussed, 20 years since the be- 
ginning of the jet/ambient interaction. The cross superimposed 
on each panel marks the best-fit temperature and emission mea- 
sure values derived from our simulated absorbed spectra. 



jet case, the EM at temperature up to a few million degrees is 
two orders of magnitude lower than in the light jet case. 

On the base of the EM(T) distribution, we expect a bright 
(I/x > 10 28 erg s _1 ) X-ray source, whose soft component (due 
to the cocoon) coul d be suppressed b y the strong interstellar 
medium absorption JFavata et alJ2002|) . We expect also that the 
X-ray emission decreases as the ambient-to-jet density ratio, v, 
decreases, leading to a brighter X-ray emission in the light jet 
case. 



3.4. X-ray emission 

From the EM(T) distributions and the MEKAL spectral code, 
we synthesized the focal plane spectra as predicted to be 
detected with the instruments on board XMM-Newton and 
Chandra (see AppendixlAlfor details), taking into account the 
interstellar absorption. To compare our numerical models with 
experimental data concerning HH 154 (the closest and best 
studied jet emitting in the X-ray band), we assumed a distance 
of 150 pc (as HH 154 that is located in the LI 551 cloud in the 
Taurus star-forming region) and an in terstellar a bsorption col- 
umn density iV H = 1.4 x 10 22 cm~ 2 JFavata et al.l2002l) . Our 
model results can be generalized to account for the other HH 
objects observations by considering different values for the dis- 
tance and the interstellar absorption. 

3.4.1 . Spatial distribution of the X-ray emission 

Assuming the jet propagating perpendicularly to the line of 
sight, from the numerical simulations we derived X-ray im- 
ages of the jet/ambient system as predicted to be observed with 
Chandra/ ACIS-I (see Appendix [Al that allows us to pinpoint 
the X-ray emission, thanks to its high spatial resolution. For all 
the models in Tab. [3] we derived the evidence that most of the 
X-ray emission produced during the jet/ambient medium inter- 
action originates from a very compact region localized at the 
head of the jet, just behind the shock front. 

Fig- ID shows an enlargement of the X-ray images as pre- 
dicted to be detected with Chandra/ ACIS-I of the head of the 
jet where most of the X-ray emission originates, 20 years since 
the beginning of the jet/ambient interaction. Note that the spa- 
tial resolution of the X-ray images in Fig. [8] is 6 times better 
than that of Chandra/ ACIS-I. In all the three cases analyzed, a 
comparison between the X-ray emitting region and the temper- 
ature and density maps in Fig.|5]shows that the X-ray source is 
coincident with the hot and dense blob discussed in Sect. 13.21 

We found that even with Chandra's spatial resolution, the 
X-ray emitting region cannot be spatially resolved (the spa- 
tial resolution of the synthesized X-ray images in Fig. [8] is 6 
times better than that of Chandra), and it will be detected as 
a point-like source. Also, significant X-ray emission is visible 
only from the hot and dense blob behind the shock front, as 
softer cocoon emission is extinguished by the the strong inter- 
stellar absorption. 

As discussed in Sect. 13.21 the X-ray emitting region for the 
three cases examined is thermally stable, and therefore the X- 
ray emission is detectable continuously during the 20-50 years 
analyzed. We found also that there are no significant varia- 
tions of the X-ray source morphology during the evolution: the 
source size varies by ±25%, always below the spatial resolu- 
tion achievable with Chandra, and there is no disappearance 
of the X-ray emitting region, according to the stability analy- 
sis. On the other hand, for some of the cases shown in Fig.|2] 
our analysis predicts a transient behaviour of the X-ray source, 
which extinguishes, few years since the beginning of the in- 
teraction between the protostellar jet and the ambient medium, 
because of the radiative cooling which dominates over the ther- 
mal conduction effects. 
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Figure 8. Synthesized X-ray emission, in logarithmic scale, as 
observed with ACIS-I, for the three cases examined, 20 years 
since the beginning of the jet/ambient interaction. At a distance 
D ks 150 pc, 100 AU correspond to about 0.7 arcsec 



From Fig.|6j the X-ray source, coincident with the hot and 
dense blob discussed in Sect. 13. 21 has a proper motion of ~ 0.7 
arcsec/yr, ~ 0.8 arcsec/yr, and ~ 2.7 arcsec/yr in the light, 
equal-density and heavy jet case respectively, assuming the jet 
axis to be perpendicular to the line of sight. In addition, we 
found that the intensity of the X-ray source decreases about 
one order of magnitude as the ambient-to-jet density contrast, 
v, decreases. Note that the heavy jet case has the higher shock 
front speed and the lower X-ray emission. 

3.4.2. Spectral analysis 

We derived the synthesized focal plane spectra as they would 
be detected with XMM-Afewfo«/EPIC-pn, characterized by an 
high effective area, with the aim to compare our model results 
with publishe d data and, in particular, with those concerning 
HH 154 dFavata et all2002l) . 



We considered two different levels of count statistics in the 
[0.3 — 10] keV band: in the low statistics case we have chosen 
an exposure time so as to obtain about 100 total photons for 
each spectrum whereas, in the high statistics case, we imposed 
about 10 4 counts for each spectrum. Albeit the latter case is 
unrealistic, given the low photon counts so far collected from 
these sources, it can help us to pinpoint some fundamental fea- 
tures of the predicted spectra. The spectral bins are grouped 
together to have at least 10 photons in the low count statistics 
case and 20 photons in the other case. 

For the models of light and equal-density jet, the synthe- 
sized spectra are well described as the emission from an op- 
tically thin plasma at a single temperature, even in the high 
count statistics case. This result is due to the strong interstel- 
lar absorption which suppresses the soft emission originating 
from the cooler plasma component in the cocoon. The best fit 
parameters derived from our simulations are shown in Tab.|4] 

The heavy jet model which best fits HH observations show 
more structured spectra than those obtained in the light jet and 
equal-density jet cases: in the high statistics case, the spectra 
are well described by a two temperature plasma emission, with 
the best fit parameters reported in Tab.|4] 

The spectral analysis reflects the structure of the EM(T) 
distribution: the spectra are sensitive to the high temperature 
portion of the EM(T) (see Fig.0, being the softer component 
suppressed by the interstellar medium absorption. On the other 
hand, in the heavy jet model, the EM(T) distribution is char- 
acterized by a bump at high temperatures broader than those 
in the other two cases (see Fig. [7J, implying that more than a 
single temperature contributes to the emission. 

Fig.|9]shows the evolution of the X-ray luminosity, Lx, in 
the [0.3 — 10] keV band, derived from the isothermal compo- 
nents fitting the spectra in the high statistics case. The light 
jet case (dots in Fig. |9} shows Lx values in good agreement 
with those observ ed in HH 154 (Lx = 3 x 10 29 erg s _1 , see 
Favat a et al . 2002). In the equal-density jet case, Lx ranges be- 
tween 7.4 x 10 28 and 2.7 x 10 29 erg s _1 , in general below the 
luminosity observed in HH 154, although its values are consis- 
tent with those detected in other HH objects (see Tab. [0- On 
the other hand, in the heavy jet case, the Lx values (crosses 
in Fig. [9} are at least one order of magnitude lower than those 
observed so far in HH objects (Tab.[0 and in HH 154 in partic- 
ular. 

3.5. Varying the jet density parameter 

As an extension of the exploration of the parameter space, we 
have varied the value of the initial jet densi t y, n;, s o far fixed to 
the value derived from Frid lund & Li seau ( 1998]) for HH 154 
(namely n- } = 500 cm -3 ). In particular, we performed the nu- 
merical simulation of a jet with initial density n- 3 = 5000 cm -3 , 
ten times denser than the ambient medium (v = 0.1), and with 
Mach number M = 20, corresponding to an initial jet veloc- 
ity Vj rs 950 km/s. We found that this model is thermally un- 
stable since the size of the X-ray emitting region at the head 
of the jet is larger than the corresponding characteristic Field 
length. As a consequence, the X-ray luminosity initially at val- 
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Table 4. Best-fit parameters for the EPIC-pn simulated X-ray spectra obtained in the low and high statistics cases, respectively, 
for the three best-fit models of Tab. [3] about 20 yr since the beginning of the jet/ambient interaction. 



Model 


counts 


7V H ± AN H 


Ti ± ATi 


EMx ± AEMi 


T 2 ± AT 2 


EM 2 ± AEM 2 


x 2 


Prob." 






(10 22 cm" 2 ) 


(10 6 K) 


(10 52 cm" 3 ) 


(10 6 K) 


(10 52 cm" 3 ) 






Light 


102 


1.4 ±0.3 


2.9 ± 1.2 


2.5 ± 10.9 






0.44 


0.88 




10317 


1.39 ±0.02 


2.9 ± 0.1 


2.2 ±0.6 






0.73 


1.00 


Equal-density 


71 


1.3 ±0.4 


3.1 ±3.8 


0.4 ±3.5 






0.49 


0.75 




10134 


1.38 ±0.02 


3.9 ±0.2 


0.33 ± 0.09 






0.72 


1.00 


Heavy 


92 


1.4 ±0.2 


6.8 ±1.7 


0.9 ±3.4 






0.39 


0.93 




9811 


1.40 ±0.04 


3.0 ±0.1 


0.2 ± 0.2 


11.8 ± 0.4 


0.016 ±0.002 


0.73 


0.99 



a Null hypothesis probability. 
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Figure 9. X-ray luminosity evolution as a function of time. 
Dots mark the light jet case, stars the equal-density jet and 
crosses the heavy jet case. The line superimposed on this figure 
mar ks Lx = 3 x 1 29 erg s _1 , the value observed in HH 154 
jFavata et al.l2002l) . 

ues w 3 x 10 29 erg/s drops of over 2 orders of magnitude in 
about 10 years. 

Higher values of the X-ray emission could be obtained in 
cases with: 1) higher initial jet density, ny, 2) higher ambient- 
to-jet density contrast, v\ 3) higher initial jet velocity, v- y 

To account for the first option (higher ny), we performed 
the same numerical simulation discussed above, but with ini- 
tial jet density 10 times greater, namely n.j = 5 x 10 4 cm -3 . 
In this case, we derived shock front velocity w 800 km/s and 
X-ray luminosity ranging between 3 and 30 x 10 29 erg/s and 
emission consistent with HH observation in general, but too 
high to reproduce the observations of HH 154. Such initial jet 
density adopted turns out to be much higher than those derived 
from observations of HH objects. IPodio et al.l 1120061) studied 
several HH objects and derived principal physical properties as 
their density. In particular they foun d that the density ranges 
between 400 and 1000 cm (see also Fridlund & Liseau 1998 
for density values in HH 154). 

As it concerns the second and third options discussed 
above, based on our exploration of the parameter space Fig. 13 
we expect that for fixed initial jet velocity, the X-ray luminosity 
increases with increasing v and vice versa. Again we conclude 
that the jet must be less dense than the ambient medium and/or 



with initial jet velocity higher than 950 km/s. Once again our 
analysis leads to the conclusion that X-ray emission originat- 
ing from protostellar jets (in particular from HH 154) is better 
reproduced by light jets with initial jet density nj = 500 cm~ 3 . 

4. Discussion and conclusions 

We presented an hydrodynamic model which describes the in- 
teraction between a supersonic protostellar jet and a homoge- 
neous ambient medium. The aim is to derive the physical pa- 
rameters of a protostellar jet which can give rise to X-ray emis- 
sion consistent with the recent ob servations o f HH objects. 

In a previous paper IBonito et alJ 120041) . we have shown 
the feasibility of the physical principle on which our model 
is based: a supersonic protostellar jet leads to X-ray emission 
from the shock, formed at the interaction front with the sur- 
rounding gas, consistent with the observations in the particular 
case of HH 154, the nearest and best studied X-ray emitting 
protostellar jet. Here we have performed an extensive explo- 
ration of a wide space of those parameters which mainly de- 
scribe the interaction of the jet/ambient system: the jet Mach 
number, M, and the ambient-to-jet density contrast, v (Fig.|2j. 
These results, therefore, provide insight in a wider set of phe- 
nomena, allowing us to study and diagnose physical properties 
of protostellar jets in general, other than that in L1551. 

One of the main results of our analysis is that only a narrow 
range of parameters can reproduce observations. The extensive 
exploration of the parameter space here discussed improves and 
extends the previous work by allowing us to constrain the main 
protostellar jets parameters in order to obtain X-ray emission, 
best fit temperature and shock front speed consistent with ex- 
perimental data. 

The range of parameters which more significantly influ- 
ences the jet/ambient evolution are shown in Tab. [2] From a 
compariso n with observ ed quantities (also shown in Tab.|2] ac- 
cording to Ballv & Reipurth 2002) it can be deduced that the 
parameters used in our model are consistent with observed val- 
ues. Note that the values of initial jet velocity are higher than 
the observed values. This apparent discrepancy is due to the 
fact that observers measure the velocity of the knots which have 
already been slowed down by the interaction with the ambient 
medium. So we need a higher initial jet velocity to account for 
these lower speed values at the working surface. However in 
the best light jet case, we derive a kinetic power L mC ch ~ 0.3 
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L , more than 2 orders of magnitude lower than the observed 
bolometric luminosity of HH 154, L^i w 40 L Q (see Tab.[Q. 
So the jet velocity values used in our simulations lead to rea- 
sonable kinetic power. Furthermore comparing the X-ray lumi- 
nosity derived in the best light jet model here discussed, we 
deduce that only a small fraction of the kinetic power is con- 
verted in X-ray emission: Lx/L me ch ~ (8 x 10~ 5 Lq)/(0.3 
Lq) = 3 x 10~ 4 . 

In the following we summarize our findings: 

- Light jet (v > 1). 

The light jet cases which reproduce the X-ray emission 
and optical proper motion observed are those with ini- 
tial Mach number M > 300 and ambient-to-jet density 
contrast 10 < v < 100. For Mach numbers lower than 
M = 300, the X-ray luminosity derived from our simula- 
tions is lower than the minimum value observed in proto- 
stellar jets ((ix)obs > 10 28 erg s _1 ) and, in some cases, 
with a transient behaviour due to thermal instabilities. The 
values M = 300 and v = 10 provide the best case which 
reproduces the HH 154 observations in terms of best fit 
temperature, emission measure and X-ray luminosity. We 
also predict a substantial proper motion. 

- Equal-density jet (y = 1). 

The three equal-density jet cases analyzed allow us to con- 
strain the initial Mach number of an equal-density proto- 
stellar jet to reproduce observations: 30 < M < 300. From 
the equal-density model with M = 100 we derive shock 
front velocity and X-ray luminosity consistent with HH ob- 
jects observations in general (see Tab.[Q. 

- Heavy jet {v < 1). 

To reproduce some characteristics of the observations, in 
the heavy jet scenario, we need initial jet Mach number 
M w 30 and initial jet density much higher than the am- 
bient density, v < 0.03, i.e. a jet 30 times denser than the 
medium or more. Lower initial Mach number in some cases 
leads to thermal instability which suppress X-ray emission 
about 5 years since the beginning of the interaction between 
the protostellar jet and the ambient medium. A jet 30 times 
denser than the ambient medium (y = 0.03) and with ini- 
tial Mach number M = 30 predicts emission from a mil- 
lion degrees plasma. Although its w s h is too high and its 
Lx is too low with respect to those observed in HH ob- 
jects in general (se e Tab. Q and in HH 154 in particular 
jFavata et alJ20 02). we cannot reject the possibility of new 
more sensitive observations which may show fainter emis- 
sion not yet detected in X-ray emitting HH objects. 

Here we discussed the best-fit models of light, heavy and 
equal-density jets in best agreement with experimental results 
from HH objects in general: a light jet with M = 300 and 
v = 10; an equal-density jet with Al = 100 and v = 1; and a 
heavy jet with M = 30 and v = 0.03. 

For each case, we analyzed the evolution of the mass den- 
sity and temperature spatial distributions derived from our 
model, the shock front proper motion and its spectral proper- 
ties, the X-ray emission and its stability. In each best-fit model 
the interaction between the supersonic protostellar jet and the 



unperturbed ambient medium leads to the formation of a hot 
and dense cocoon surrounding the jet and smoothed by the 
thermal conduction. Just behind the shock front, there is a hot 
and dense blob from which the harder and bright X-ray emis- 
sion originates: the strong interstellar absorption suppresses the 
softer component due to the cocoon. In all cases examined, the 
X-ray emitting region is thermally stable, i.e. thermal conduc- 
tion effects prevent the collapse of the source due to radiative 
cooling, and shows a detectable proper motion. 

To compare our findings with HH 154 observations, we 
have rejected any equal-density and heavy jet case which 
show too high shock front velocity (v s h, ~ 600 km s _1 and 
w 1900 km s~\ respectively) and too low X-ray luminosity 
(Lx ~ 10 29 erg s _1 and « 10 28 erg s~ x ) with respect to v s h 
and Lx values observed in HH 154 (w s h ~ 500 km s _1 and 
L x = 3 x 10 29 ergs" 1 ). 

In the best light jet case, we derived a particle density n ~ 
6500 cm~ 3 and a velocity v s h ~ 500 km s _1 for the X-ray 
emitting region, leading to a momentum mv = 2 x 10 -6 M© 
km s _1 . This value is consistent with the u pper limit 7 x 10~ 4 
A1 Q km s" 1 obtained for HH 154 dFridlund & Liseaulll998l) . 
This leads to the conclusion that the protostellar jet cannot drive 
the molecular outflow, whose momentum has been e stimated 
to be between 0.15 and 1.5 A/ Q km s _1 ( Fri dlund & Liseaul 
1998). The relations for the mass loss rate and the mechanical 
luminosity 



M = 1.8 x 10" 9 M Q /yr 
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( 3 ) 

V 500cm" 3 / 



( V J 
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5 x 10 14 cm / V 500cm- 3 / V 1400fcm / 



lead, respectively, to values 3 an d 2 orders of magnitude lower 
than expected in CO outflows JCabrit & Bertoutlll992l). This 
resul t supports the conclusion discussed by Fridlund & Liseau 
( 1998) that the jet origin is probably different from that of the 
CO outflow. 

We estimated the values of the momentum, mass loss rate 
and mechanical luminosity for each model which better re- 
produce the X-ray observation of HH objects shown in Tab. 
We derived: M - 10~ 9 Af /yr; L moch = (0.3 - 1.7) 
L Q ; mv < 10 -6 M km s _1 . These values are several 
orders of magnitude lower than those ob served in CO out- 
flows (in HH 2, HH 154 and HH 80/8 1 . iMoro-Martm et all 
Il999t iFridlund & Liseaulll998t lYamashita et al.lll989l) . From 
the comparison between our results and the observations of CO 
outflows, we conclude that the simulated protostellar jets which 
best reproduce X-ray observations cannot drive molecular out- 
flow. 

On the base of our analysis, we conclude that the light jet 
scenario, 10 times less dense than the ambient medium (y = 
10) and with initial Mach number M = 300, with v sh w 500 
km s _1 and Lx ~ 3x 10 29 erg s _1 , is the best case to reproduce 
the HH 154 observations. This is also supported by the optical 
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Figure 10. Mean value of the X-ray luminosity as a function 
of the mean value of the best fit temperature, derived from our 
model synthesized spectra. Crosses correspond to Mach num- 
ber M = 1000, stars to M = 300, diamonds to M = 100, 
triangles to M = 30 and squares to M = 10. Increasing size 
of each symbol corresponds to increasing value of the ambient- 
to-jet density contrast, v = 0.01-300. The shaded region refers 
to the (Lx, T) values consistent with observations of HH ob- 
jects in general. 

observations of HH 154 jFridlund et alJl2005l) from which the 
light jet scenario can be deduced, according to the Hartigan 
( 1989) model. 

More in general, in Fig.^]we show the values of Lx vs. T 
derived from the spectra synthesized from our model as a func- 
tion of different values of M and v. Crosses mark cases with 
M = 1000, stars with M = 300, diamonds with M = 100, 
triangles with M = 30 and squares with Al = 10. Bigger 
symbol sizes correspond to higher values of the ambient-to-jet 
density contrast, v, in the range 0.01 to 300, as in Fig. 13 The 
shaded zone marks the range of parameters consistent with ob- 
servation (Tab.Q. We have chosen Lx one order of magnitude 
lower than the minimum observed so far to account for fainter 
sources. 

For a fixed value of M, Lx and T decrease with increasing 
v and, for a fixed value of v, Lx and T increase with increasing 
Al. From the figure it is possible to derive Al, v (velocity and 
density) of the protostellar jet to be compared with observations 
(in terms of Lx and best fit T). Predictions about fainter not yet 
discovered sources can also be made. 

Furthermore the results derived from the variation of the 
initial jet density parameter, nj, discussed in Sect. 13.51 lead to 
the conclusion that, even with different initial density values, 
a protostellar jet must be less dense that the ambient medium 
and with high initial velocity (1000 km/s or more) in order to 
reproduce HH objects observations. 

Our model predicts in all cases a significant proper motion 
of the X-ray source, with values which, in the case of HH 154 
would be measurable with Chandra, providing a clear test of 
the model scenario. 

As discussed bv lFavata et alJ J2006h . a 100 ks observation 
was performed in 2005, showing, when compared with the 
2001 observation, a more complex scenario, i.e. both a mov- 



ing and a stationary source were detected in HH 154, giving the 
source, in 2005, a "knotty" appearance. Thus, while a traveling 
shock (likely based on the basic physics explored in the present 
work) is apparently present in HH 154, the source structure is 
more complex. 

The comparison between our model of a continuous super- 
sonic jet through an unperturbed surroundin g med ium and the 
new Chandra data, discussed in lFavata et alJd2006h . shows that 
the model reproduces most of the physical properties observed 
in the X-ray emission of the protostellar jet (temperature, emis- 
sion measure, etc.). At the same time, it fails to explain the 
complex evolving observed morphology, showing, most likely, 
that the jet is not continuous. 

A possible scenario, which we will test in the future, is 
based on similar physics, but in the presence of pulsating jets 
(instead of the continuous jet here examined), or alternatively 
on interaction between the jet and an inhomogeneous ambient 
medium which can lead to the knotty structure observed inside 
the jet itself. We also plan to explore other physical mecha- 
nisms, different from the moving shock at the tip of a super- 
sonic jet, stimulated by the above mentioned Chandra observa- 
tions of HH 154, such as steady shocks formed at the mouth of 
a de Laval nozzle. New observations of the evolution of HH 154 
will however be necessary to understand the phenomenon and 
to further constrain the model scenario. 

The continuous jet model discussed here is a useful and 
necessary building block toward more complex (e.g. with a dis- 
continuous time profile) models. The X-ray emission from a 
pulsed jet, for example, will still take place at the shock front, 
and will thus be based on the same physical effects and prin- 
ciples as observed in a continuous jet. Our (simpler) continu- 
ous jet model is a very useful tool to infer the right parameter 
values to use in future, more complex models needed to also 
reproduce the observed morphology: using our continuous jet 
model results, shown in Fig. [2] and in Fig. it is possible 
to derive the initial jet velocity and ambient-to-jet density ratio 
needed for the new bullets to produce to X-ray luminosity and 
best fit temperature consistent with observations. 

A limitation of our hydrodynamic model is the hypothesis 
of pressure balance between the jet and the ambient medium, 
in order to obtain the observed jet collimation. Most models 
of jet collimation suggest the presence of an organized ambient 
magnetic field which is known to be effective in collimating the 
plasma. As a follow-up of our analysis, we are developing an 
MHD model of protostellar jet that will allow us to relax the as- 
sumption of an initial pressure equilibrium. The comparison of 
the MHD model results with the X-ray observations will pro- 
vide a fundamental tool to investigate the role of the magnetic 
field on the protostellar jet dynamics and emission. 
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Appendix A: Synthesizing the X-ray spectra 

From our 2-D numerical simulations, we synthesized the ab- 
sorbed focal plane spectra to be compared with observations 
by using the following procedure. 

As a first step, from the integration of the hydrodynamic 
equations 0and E| we derive the temperature and density 
2-D distributions in the computational domain. 

We reconstruct the 3-D spatial distribution of these phys- 
ical quantities by rotating the 2-D slabs around the symme- 
try axis. Then we derive the emission measure, defined as 
EM — J n e rindV (where n e and nn are the electron and hy- 
drogen densities, respectively, and V is the volume of emitting 
plasma). 

From the 3-D spatial distributions of T and EM, we derive 
the distribution of emission measure EM(T) for the compu- 
tational domain as a whole or for part of it: we consider the 
temperature range [10 3 — 10 8 ] K, divided into 74 bins equis- 
paced in log T; the total EM in each temperature bin is ob- 
tained summing the emission measure of all the fluid elements 
corresponding to the same bin. 

From the EM (T), using the MEKAL spectral code 
jMewe et alJll985l) for optically thin plasma, we derive the 
number of photons in the i-th energy bin as 



where D is the distance of the object, Ei is the energy in the 
i-th bin, P(Tk, E) describes the radiative losses as a function 
of the energy and of the temperature in the k-th bin. 

To compare our model results with observations, we syn- 
thesize the focal plane spectrum, Ci, as predicted to be ob- 
served with the Chandra/ ACIS-I or XMM-Afewfew/EPIC-pn X- 
ray imaging spectrometers taking into account the spectral in- 
strumental response: 



where t exp is the exposure time, A(E) is the effective area and 
M(i, E) is the instrumental response. 

Finally we take into acc ount the interstellar me dium ab- 
sorption column density, JVh JMorrison & McCammo nl 19831) . 
and we analyze the absorbed focal plane spectrum with XSPEC 
VI 1.2 in order to compare our findings with published experi- 
mental results. 




(A.l) 




